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We derive analytical expressions for the correlation functions of the electronic conductance fluctuations of an 
open quantum dot under several conditions. Both the variation of energy and that of an external parameter such 
as an applied perpendicular or parallel magnetic fields are considered in the general case of partial openness. 

These expressions are then used to obtain the ensemble averaged density of maxima, a measure recently sug¬ 
gested to contain invaluable information concerning the chaoticity of the system. The correlation width is then 
calculated for the case of energy variation and a significant deviation from the Weisskopf estimate is found in 
the case of two terminals. The results are extended to more than two terminals. All our results are analytical.The 
use of these results in other fields, such as nuclei, where the system can only be studied through a variation of 
the energy, is then discussed. 

PACS numbers: 05.45.Yv, 03.75.Lm, 42.65.Tg 


INTRODUCTION 

Chaotic quantum systems, such as open Quantum Dots fD, 
graphene flakes 03, and nuclei exhibit common universal 
features characterized by fluctuating observables, such as the 
electronic conductance in the former two, and the compound 
nucleus cross section in the latter SO- Useful information 
about the system in this regime of the dynamics are obtained 
through the average of the observables and through the cor¬ 
relation functions defined as the average of products of two 
fluctuating observables at two values of the independent vari¬ 
able, such as the energy or an external field. The correlation 
function is of paramount importance as it measures the degree 
of coherence present in the otherwise fully chaotic system. 
This measure resides in the correlation width which specifies 
the shape of the correlation function. In the case of paramet¬ 
ric energy variation the shape is Lorentzian, as shown more 
than half a century ago by Eric son m, while in the case of a 
varying external parameter such as a magnetic field, the shape 
is a square Lorentzian, EHID. In CD, the effect of tem¬ 
perature on the conductance correlation function was inves¬ 
tigated. Deviations from these shapes have also been inves¬ 
tigated OH E). For their importance in the study of chaotic 
mesoscopic systems, nanosystems, resonators and in nuclei, 
it is clear that a detailed investigation of the correlation func¬ 
tions and their potential deviations from the above mentioned 
two extreme shapes and the resulting correlation widths is 
called for. In this paper we calculate analytically the corre¬ 
lation functions for the general case of partially open system. 
Three types of external magnetic fields are considered. The 
energy correlation function is discussed in details and the cor¬ 
responding correlation width is extracted and found to be sig¬ 
nificantly different from the Weisskopf estimate. This last re¬ 


sult is quite important as it affects the value of the dwell time 
in these systems. 

QUANTUM CHAOTIC SCATTERING 

Quantum Chaotic Scattering (QCS) is a widely occurring 
phenomenon in physics. It operates in a variety of systems, 
such as, electronics m, spintronics 03, biomolecules 
m, disordered mesoscopic nanostructutes (33, and the 
compound nucleus ffiD. The emergence of the phenomenon 
is, just as in other cases of quantum chaos, directly related to 
the intrinsic chaotic dynamics associated with quasi-bound 
states of a quantum system. Random Matrix Theory (RMT) 
supplies the formal S-matrix of QCS, as it describes the 
intrinsic Hamiltonian which governs the dynamics and the 
scattering. The empirical manifestation of QCS is evidenced 
by universal fluctuations observed in electronic conductance 
in open Quantum Dots, and in Graphene flakes, transmittance 
in microwave and acoustic resonators, and in compound 
nucleus cross sections. 

The S-matrix describing QCS is given by 

S(s,X) = t~2mW^(£-H(X) + iTrWW^)- 1 W , (1) 

where H(X) is a random Hamiltonian matrix of dimension 
M x M that describes the resonant states in the chaotic 
system, which is subject to the influence of an external 
parameter, X. The number of resonances is very large 
(M -A oo). The matrix W of dimension M x (N) contains 
the channel-resonance coupling matrix elements. Using the 
above S-matrix, one is then able to calculate observables. 
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The ensemble average (5 CC '(£)5/ C , (e')) supplies the average 
conductance or cross section (e = s'), while the four-point 
function, (5 cc /(£i)5/ c ,(£i)5 cc '(£ 2 )5/ c ,(£ 2 )), furnishes the 
correlation function, characterized by the correlation width, 
of great importance in the study of chaotic quantum systems. 

In application to conductance statistics in Quantum Dots 
with three or more terminals it is convenient to represent the 
5-matrix as, 


5 = 


r 11 

^12 

h 3 

t' 

t 12 

T22 

t23 

f< 
c 13 

f' 

c 23 

f33 


( 2 ) 


where tij indicates the probability amplitude of transmission 
of the channel(s) contained in the terminal i for the channel(s) 
contained in the terminal j and r n denotes the probability am¬ 
plitude of reflexions of the channels in the terminal i. 

The correlation function obtained from the above is invari¬ 
ably of a Lorentzian shape, in the case of a variation of the 
energy, or a square Lorentzian, in the case of a variation of the 
external parameter, X. Deviations from these limiting cases 
are expected and can be derived ED El- Calculation of the 
ensemble averages, is rather difficult, and only the case of the 
Gaussian Unitary Ensemble average has been performed an¬ 
alytically. Most applications of the above 5-matrix, involves 
numerical simulations using random matrix generator. To ob¬ 
tain analytical results, an alternative method has been devised, 
based on the distribution of the 5-matrix itself, which, being 
unitary, follows the Dyson circular ensemble. The Stub model 
is such an alternative EHE]. It involves the use of the fol¬ 
lowing parametrized form of the 5-matrix. 


we are assuming the equivalent coupling for channels 
in the same terminal. The 77-matrix (representing the 
external fields) is the stub model counterpart of order 
(M — N x M — N) as described by (T9). The parameter X 
represents, e.g., the type of external applied magnetic field 
employed. The quantity T = fr(fi 2 ij 2 ). 

The calculation of the ensemble averages can be done in a 
closed form using the stub 5-matrix above. Details of such 
calculations can be found in the original references. The en¬ 
semble average of the conductance, Q = tr(t 12 ^ 12 )’ is easily 
calculated, lf20l . following the Ref. fl9l of a chaotic QS. We 
obtain 


(S) 


T 1 T 2 [ 2-/3TiT 2 +T 2 Ti 

T [ p T 2 


(4) 


where T, = NiTi is the total transmission coefficient for the 
terminal i with equivalent channels and T = Ti + I 2 + 75, 
and the e-dependence disappear in the ensemble average. In 
the case of the compound nucleus, the quantity Q becomes 
the cross section and the its ensemble average , similar to eq. 
(4), is known as the Hauser-Feshbach cross section ED, with 


the factor 


1 - 


2-/9 T1T2+T2T1 

/3 T 2 


representing what is known 


as the elastic enhancement factor, which is known exactly at 
arbitrary overlap both for orthogonal and unitary symmetry, 
E 2 E 1 . Note that the number of channels, N, in all the 
terminals is taken to be large, N » 1, and thus we are in 
the semiclassical regime. Note further that the symmetry pa¬ 
rameter /3 takes on the values 1, in the Gaussian Orthogonal 
Ensemble, 2 in the Gaussian Unitary Ensemble, and 4 in the 
Gaussian Symplectic Ensemble. 


THE CORRELATION FUNCTIONS 


THE STUB MODEL 

We assume that the particles dynamics is ballistic and er- 
godic, and we model the system statistical properties using 
the random matrix theory (RMT)-based stub model. Follow¬ 
ing references Emu, for particles with spin, the scattering 
matrix S can be represented as a unitary matrix with quater- 
nionic entries and is parameterized as 

S(e, B) = S + VU[1 - /C t 77(£, (3) 

Here, 77-matrix, M x M, is the scattering matrix counterpart 
of an isolated quantum system, while S is the average of 
the scattering matrix of the system 5, which has dimension 
N x TV. The M stands for the number of resonances of the 
system, while N = Ni + N 2 + -/V 3 is the total number of 
open channels. The universal regime requires M 2 > N. The 
KL- matrix is a projection operator of order (M — N ) x M , 
while V, of order N x M, describes the channels-resonances 
couplings. Their explicit forms read JCij = <5j _'Pij = 
di ag (( fi’Lj 'Ji 1 , ;V| .j Ni -\- N ‘2 .j \/7)) and 5/ j 

diag(<5i,j-\/l — 7\, j-\A — T 2 , Si + N 1 +N 2 ,jV^ — T 3 ), 


As for the correlation functions, the calculation using the 
stub model can be carried out as done in l24j . In the following 
we merely write down the correlation functions for a chaotic 
Quantum Dot in the case of only two terminals, specified by 
7’| and T 2 and and consider Ti = T 2 = T, we find for the 
energy variation we find, 

C(SE) 3T(2 - T) - 2 4 [1 + T(T - 2)] 

1 / 8/3 ~ 1 + (SE/T e Y + [1 + {SE/TeY) 2 

where /3 is a measure of the universality class to which H 
pertains. It should be mentioned that magneto-conductance 
correlation functions, were calculated using the Hamiltonian 
approach, Eq. (1), and the study of a relevant asymptotic ex¬ 
pansion in inverse channel number was reported in l25l . 

For the variation of the external parameter, we consider 
three cases. An external perpendicular magnetic field, B ±, 
and external perpendicular magnetic filed acting on the spin- 
orbit interaction, the so-called Rashba-Drasselhaus field, II± 
EBE3, and a parallel magnetic field, 7//, whose effect has 
been studied recently by ll28l . 
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FIG. 1; Fig.(la) show a typical Lorentzian behaviour of the correlation function and the correlation length T corr associated with a typical transport observable between the terminals 
1 and 2 probed by another non-linear terminals 3, as the inset indicates. Fig.(lb) shows a deviation from the Weisskopf length in the semiclassical regime occasioned by the presence 
of another terminals. Fig. 1(c) shows the ratio D = r corr /Tw as a function of T 3 exhibiting the attainment of the Weisskopf length only for a completely open (T 3 = 1) or 
completely closed (T 3 = 0) extra-terminals lumped into an effective one. See text for details.. The continuous line is the analytical result and the dots are the results of the numerical 
Hamiltonian simulation using the S'-matrix of Eq. (1). Fig.(ld) shows the deviation of the Weisskopf length analytically and numerically. 


[eq.(14) of the present work] were analyzed in great detail 
C{6X) 2T(1-T) 2 + T(3T — 4) there as well. 

1/8/3 ~ 1 + (6X/T x ) 2 + [l + {SX/TxY} 2 


3 T 2 4- 4 — 4 T T 2 — 2 T 

Finally for the case of a parallel magnetic field considered re¬ 
cently in lf28l 

Cm\ |) = T(2~T)~l 
1/8/3 l + iSH^/T^) 2 

This last correlation function is new and deserves some 
discussion. It does not depend on the openness parameter, 
the transmission coefficient, T, and it has a pure Lorentzian 
shape in contrast to the cares of perpendicular magnetic field. 

It should be mentioned that a thorough study of the energy 
dependence and correlation length in quantum chaotic scat¬ 
tering with many channels was performed in lf32l A general 
relation was established there between fluctuations in scatter¬ 
ing and the distribution of complex energies (poles of the S 
matrix). In particular, the correlation length was shown to be 
given by the spectral gap in the pole distribution, and the devi¬ 
ations of the gap from the (semiclassical) Weisskopf estimate 


THE AVERAGE DENSITY OF MAXIMA 


The important feature that characterizes these correlation 
functions is that they all (except the case of H\\) deviate from 
pure Lorentzian or square Lorentzian shape, for an arbitrary 
value of the openness probability. The application to the case 
of the compound nucleus, Eq. ( 5), allows the investigation 
of the statistics of resonances both in the weak (isolated reso¬ 
nances) and strong (overlapping resonances) absorption cases, 
as well as in the intermediate cases. Several quantities can be 
obtained from the correlation functions. The average density 
of maxima in the fluctuating observable, is one of them. In 
Ref. |f30l , this quantity was derived and analyzed for C\E) 
and C{X). For completeness we give below the main results 
and extend them to other types of applied magnetic fields. 
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For the energy variation. 


(9) 


7e) 


73 19 T 2 - 18T+ 10 

ttTb V 5T 2 - 10 T + 6 


73 


( 10 ) 


For the case of an external perpendicular magnetic field. 


, x 73 / 7 T 2 - 10T + 6 

{PX ' ~ T27rr A - V 2T 2 -3T + 2 

For the Rashba-Dresselhaus field, 


( 11 ) 


73 )7T 2 - 10T + 6 

2nT h± V 2 T 2 - 3 T + 2 


( 12 ) 


Interestingly for the case of a parallel magnetic field, the re¬ 
sult is independent of the openness parameter and is identical 
to the Brink-Stephen m one, 


{PH\\ ) 


73 

nT Hll 


0.55 


(13) 


values that reach up to 1.1. In our work here, we can read out 
the change in the correlation width as, 


n _ T ^rr _ / 5T 2 — 10T + 6 

T w V 9T 2 - 18T+ 10 ‘ 

For almost closed system, T « 1, the deviation reaches the 
value D = 73/5 = 0.81. Of course in the other limit of a 
completely open Quantum Dot, T « 1, the ratio I) attains the 
value of unity as expected. . 

It is interesting to generalize our result to the case of two 
terminals in the presence of other terminals lumped together 
as an effective terminal characterized by T 3 , and study the 
variation of the correlation width as a function of X 3 . Thus, 
we extend the stub method for the new calculation of the 
general correlation function, including the quantum interfer¬ 
ence terms. Using the same stub 5-matrix described above, 
together with the typical large number of diagrams for the 
correlation function 03 , we calculate the averages of prod¬ 
ucts of two observables, Q(E)Q(E'), and obtain after setting 
71 = T 2 = N, and N 3 = N 2 = N 3 = N. Note that in this 
case of more than two terminals, the result we have obtained 
are valid for T) = T 2 = 1, to get. 


C(SE) 

I//? 


and, 


[A(T 3 )-87V 2 7?(1-T 3 ) 2 ] 
N 4 (2 + T 3 ) 6 [1 + (5E/T W ) 2 } 

, T 3 A 2 (T 3 ) 

2/3 N 3 (2 + 7 3 ) 6 [1 + (5E/T W ) 2 } 

8T 3 2 (i - 2 r 3 + rf) 

(2 + T 3 ) e [l + (SE/T w ) 2 } 2 


THE CORRELATION LENGTH AND WEISSKOPL’S 
ESTIMATE 

The important point to emphasize here is that both the cor¬ 
relation function and the average density of maxima are char¬ 
acterized, for a fixed value of the openness probability, by a 
single quantity, the correlation width. In the energy variation 
case, this width is the inverse of the dwell time and is usually 
estimated using the Weisskopf expression lf33l . 

t e *t w = ^Y,t c ( 14 > 

c 

where, A is the average spacing between the resonances in the 
chaotic system, and the sum extends over all the open chan¬ 
nels reached through the transmission coefficient, T c . Devia¬ 
tion from the Weisskopf estimate was calculated in ll34l using 
the 5-matrix of Eq. (1). With the help of a random matrix 
generator, these authors calculated numerically the transmit¬ 
tance correlation function for microwave resonators, and ob¬ 
tained the correlation width, T corr as the width at half maxi¬ 
mum. They found for the ratio D = T corr /T-w vs. ^ C T C , 


-4r(r 3 ) 

iV 4 


A 2 {T 3 ) 

N 3 


Ai(T a ) = T 4 + 2(2 p + T 3 )( 2 + T 3 )T 3 3 

[2Tg - 4(2 + T 3 ) 2 T 3 

16(2 + T 3 ) 2 ] T 2 

2(2 p + T 3 ) [—4T 2 (2 + T 3 ) 

12(2 + T 3 ) 2 ] T 3 + 8(2 + T 3 ) 4 

A 2 {T 3 ) = T 3 (2 + T 3 ) 2 [2(2 + T 3 ) + 1] 

(17) 


The Eq.(p~ 6 ]> shows convincingly that the correlation func¬ 
tion is not a Lorentzian on the Weisskopf width scale, which 
is violated and no longer the transport correlation width. The 
Fig. (lb) shows the modification of the Weisskopf length as 
a function of the T 3 . On the another hand, the amplitude of 
the universal fluctuations (variance) var[Q\ = C(0) = 1/(3 is 
maintained regardless of the additional terminals. The Eq.([T 6 |) 
can be written in Lorentzian form, using some algebra, as in 
the following 


C(6E) 


1 


1 + {5E/T corr ) 


2 ’ 


r corr — r \y D 


(18) 
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with 


D 


\J 64Tg (1 - T 3 ) 4 + [A^Tz) + 5 20 T 3 A 2 (T 3 )} 2 


Ai(T 3 ) + S 20 T 3 A 2 (T 3 ) 

8 T 3 (1 — T 3) 2 
Ai(T 3 ) + <5 2 ^T3A 2 (T 3 ) 


(19) 


which is highly non-liner as a function of T 3 . The form of the 
Eqs.([T 8 | and ( fl9| ) show analytically the effect of the presence 
of more terminals on the correlation properties of the other 
two terminals. It simulates absorption of the flux in the two- 
terminal subsystem. The non-linear effect disappears in the 
limits T 3 = 1 (ideal) and T 3 —» 0 (opaque) for which D —> 1, 
and the correlation length approaches the Weisskopf length, 
as Fig.(lc) shows clearly. These results were verified by a 
numerical simulations with random matrix generator using the 
S'-matrix of Eq. (1). 
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FIG. 2: Typical transport observable Q as a function of e/Yw for 
N = 50 open channels coupled to 5 ■ 10 3 resonances. For a sin¬ 
gle realization of H in both fluctuations and the same observable 
Q, the top indicates a single measure in the absence of the extras- 
terminals and the bottom indicates extra-terminals with non-linear 
coupling T 3 = 0.3. 


We perform a numerical simulation using the Hamiltonian 
model of Eq. (1) with a configuration of 50 open channels and 
5 • 10 3 resonances in an energy range Ae /7 £ [—100,100]. 
As shown in Fig. [2], the transport observable Q is “experi¬ 
mentally” obtained for the same scattering QS. For such sin¬ 
gle QS the Q is largely affected in the presence of the extra 
non-linear terminals in two forms. Firstly, the reference line 
of fluctuations is suppressed (from the top to the bottom plot 
of Fig.[2]), as can be expected using Eq.(4). Secondly, the 
number of maxima increases from 99 to 101 in the same en¬ 
ergy interval, confirming nicely our analytical findings. A sin¬ 
gle experimental trace can confirm the violation of Weisskopf 
width in the presence of the extra terminals. 

Thus, we find excellent agreement between our analytical 
calculation and the numerical one. In a way we are also con¬ 
firming the results of l34l obtained for transmittance in a mi¬ 
crowave cavity. 


CONCLUSIONS 

In this paper, we analytically calculate the correlation func¬ 
tions for chaotic systems using the Quantum Chaotic Scatter¬ 
ing theory. We show that in the case of energy variation and 
in the variation of an external magnetic field, these functions 
deviate significantly from the expected Lorentzian and square 
Lorentzian shapes. The parameter that measures this devi¬ 
ation is identified as the transmission coefficient which acts 
as the ’’openness” probability. We further identify the devi¬ 
ation of the correlation width from the Weisskopf width for 
an arbitrary observable of the quantum transport. The stub 
calculation of the ensemble average of the products of four 
scattering matrices allows finding this deviation, of the or¬ 
der of 10 %, and reveals the topological effects of other ter¬ 
minals on the chaotic quantum transport. The consequences 
of our results include the increase of the dwell time of an 
arbitrary chaotic scattering system. The results are general 
and applicable to any scattering amplitudes between two ter¬ 
minals within the diagrammatic approach of the stub model, 
which is an 1/N expansion method. Therefore the deviation of 
the correlation width from the Weisskopf estimate can occur 
in spin and/or charge channels for electronic nanostructures, 
the transmittance of antennas, sub-lattices and/or sub-valleys 
channels for graphene flakes li35l . etc. Technological applica¬ 
tions of the results includes artificial atomic clocks by means 
of dwell time modifications and cryptography on transport ob¬ 
servables. 
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